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We consider a time dependent trap externally manipulated in such a way that one of its bound 
states is brought up towards the continuum threshold, and then down again. We evaluate the 
probability P stay for a particle, initially in a bound state of the trap, to continue in it at the end 
of the passage. We use the Sturmian representation, whereby the problem is reduced to evaluating 
the reflecting coefficient of an absorbing potential. In the slow passage limit, P stay goes to 1 for a 
state turning before reaching the continuum threshold, and vanishes if the bound state crosses into 
the continuum. For a slowly moving state just "touching” the threshold P atay tends to a universal 
value of about 38%, for a broad class of potentials. In the rapid passage limit, P stay depends on 
the choice of the potential. Various types of trapping potentials are considered, with an analytical 
solution obtained in the special case of a zero-range well. 

PACS numbers: 03.65.-w 


I. INTRODUCTION 

Recent technological developments have renewed the 
interest in the dynamics of a particle, or particles, 
trapped in bound states of time dependent potentials. 
External manipulation of Hamiltonians with both dis¬ 
crete and continuum spectra routinely occur in applica¬ 
tions such as metrology and quantum information pro¬ 
cessing. The presence of a continuum plays an important 
role in atom lasers mm, in the preparation of atomic 
pulses with a known velocity distribution |3|, or in the 
production of few-body number states HHEj. Quite often 
a continuum is responsible for undesirable loss of trapped 
particles, as it happens in transport of trapped ions, or in 
trapped ion atomic clocks. An obvious way to avoid such 
loss is to manipulate the trapping potential sufficiently 
slowly (adabatically) so that the trapped particle would 
remain trapped throughout the evolution. 

The question of adiabaticity in bound-to-continuum 
transitions, studied by various authors, (SHE] leaves 
room for further discussion, even in regard to its for¬ 
mulation. As a trapping potential becomes shallower, 
a bound states is brought closer to the continuum, and 
eventually joins it. With this drastic reorganisation of the 
adiabatic spectrum, application of methods developed for 
levels crossing situation, such as the original Landau- 
Zener model m and its numerous generalisations, is at 
best problematic. Moreover, in an experimental situa¬ 
tion one is likely to control the shape of the trap, so that 
the evolution of the energy of the bound state near the 
continuum threshold, must be deduced from that of the 
potential. With this in mind, one may be interested in 
asking two distinct questions. Firstly, let the depth of 
the trap decrease linearly in time. When the evolution 
stops, what is the probability to remain in the modified 
bound state? Secondly, let the depth of the potential first 
decrease, and then increase again, e.g., being a quadratic 


function of time. What is the probability to remain in 
a bound state at the end of the passage? The first case 
was studied in m. In this paper, we consider the second 
generic case, where a time dependent trap is manipulated 
in such manner, that a bound state completes a passage 
near the continuum threshold, first rising towards it, and 
then moving away again. There are three possibilities: 
the state may "turn” and begin the downward part of 
its journey before reaching the threshold. Alternatively, 
it can just ’’touch” the threshold once, or cross into the 
continuum temporarily, to reappear at a later time. In all 
cases we will want to know the probability for remaining 
in the initial state, or, more generlally, inside the well, 
once the passage is completed. 

As in |14J we will employ the Sturmian technique, de¬ 
veloped in Refs. fTobTSl for applications in the theory 
of atomic collisions. In this way, we reduce the prob¬ 
lem of solving a time-dependent Schroedinger equation, 
to a simpler problem of determining the reflection coef¬ 
ficient of a complex valued ’’potential”. This, in turn, 
will allow us gain further insight into what happens near 
a continuum threshold, and occasionally obtain an exact 
analytical solution to the problem. 

The rest of the paper is organised as follows: is Sect. 
II we will formulate the problem of a time dependent 
trap, which can lose a previously bound particle to the 
continuum. In Sect. Ill we introduce the Sturmian basis, 
and use it to expand the particle’s state. In Sect. VI we 
consider a zero-range well, and formulate the adiabatic 
condition for the passage. In Sect. V we solve the zero- 
range problem exactly for the case where the bound state 
just touches the continuum threshold. We will show that 
the probability to remain in the well is independent of 
the rate of change of the potential, and always equals 
approximately 38%. In Sect. VI the general case for a 
zero-range potential is analysed. In Sect. VII we consider 
the Sturmian representation for a rectangular potential, 
and the corresponding adiabatic limit. In Sect. VIII we 
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employ the single-Sturmian approximation in order to 
describe the particle’s evolution in a rectangular well. In 
Sect. IX we show that the 38% rule of Sect. V applies 
universally in the slow passage limit to a wide class of 
potentials whose evolution is quadratic in time. Sect. X 
contains our conclusions. 


II. LOSS AND RECAPTURE OF PARTICLES 
BY A TIME DEPENDENT POTENTIAL WELL 

We start by considering a particle of a mass p in a 
one-dimensional potential, 

K 

V(x,t) = -W(x)Y J V ik h k , (1) 

fc =0 

where W(x) is normalised by the condition 
f° c xj W(x)dx = 1. The potential is obtained by 
varying the magnitude of a finite-range potential well 
—W{x) < 0 by means of a time dependent factor, so that 
whenever Y^k= 0 V^ k H k turns negative, V(x,t) becomes 
a barrier, which doesn’t support bound states. The 
question we will ask is the following one: if a particle is 
put into one of the bound states of the well, <j> n , what is 
the probability to still find it there at some time in the 
future? The Schroedinger equation (SE) to be solved 
has the form (we will use Ti = 1) 

id t ^>(x,t) = -dl^f/2/j + V(x,t)^>, (2) 

and we will assume that the potential is a deep well in 
the distant past and future, V(x,t) < 0 for t —> ±oo. A 
particle in a bound state <f> m (x % t), = 1 with a 

large negative energy E m (t) < 0 should continue in it for 
some time, before approaching the continuum threshold 
mi- For T(a;,t) in Eq. ([3]) we, therefore, write 

lim ^f(x,t) = exp[-i [ E m (t')dt']<j> n (x,t), (3) 

t—> — OO J 

Similarly, for t —► oo, we should have 

lim ^>(x,t) = (4) 

t->-oo 

y A mn exp[-7 f E n {t')dt'](j) n {x,t) + 5^f(x,t), 

n 

where the first term corresponds to the particles which 
remained in the well, although possibly not in the same 
state, and (54/ describes the particles lost to the contin¬ 
uum during the passage. Thus, the total probability for 
the particle to stay in the well is given by 

p% av = Y, p ™=Y, i a -i 2 . ( 5 ) 

n n 

In the following will will consider the simplest case of a 
passage, which is quadratic in time, 

V(x,t) = {£ — v 2 t 2 )W(x), (6) 


continuum 



t l 

t < 0 t > 0 


FIG. 1. (Color online) Schematic diagram showing the evolu¬ 
tion of the potential well described by Eq.|6]|. At first the well 
is becoming shallower, thus bringing its bound state (thick 
solid) closer to the continuum. Later the well deepens, bring¬ 
ing the state down, and possibly bringing more bound states 
into the well. At t = 0, the bound state may still exist, if 
£ < 0, ’’touch” the continuum threshold, if £ = 0, a disap¬ 
pear for a while, if E > 0. 

and of a particle trapped in an ascending bound state 
in the distant past, which may remain trapped in one of 
the descending states, or be ejected into the continuum as 
t —> oo (see Fig.l). In particular, the ground state of the 
well (which in one dimension exists as long V(x, t) < 0 
[13]) will ’’turn” before reaching the continuum threshold 
if £ <0, ’’just touch” it if £ = 0, or, for £ >0, disappear 
at t = —VE/v, before re-appearing again at t = \/£/v. 
In all three cases, we will be interested in the probabilities 
Pmn V defined in Eq. (JsJ). 

III. STURMIAN EXPANSION OF THE TIME 
DEPENDENT STATE 

With the help of the Fourier transform, 

^(x,t) = J dwexp(—*wt)4/(a;, w), (7) 

we rewrite Eq. © as 

K 

uj^{x,oj) = -dl^/2ii-W{x)Y^(-i) k V {k) dt i &. ( 8 ) 

k =0 

and look for a suitable basis in which to expand ^(Xjw). 
Using the set of the positive-energy scattering states de¬ 
scribing particles incident on V (x, t) from left and right is 
one option, yet there is a more convenient one. Particles 
ejected from the well should be described by outgoing 
waves on both sides of the potential. Sturmian basis sets 
with the desired properties are well known in literature 
m They are obtained by imposing outgoing boundary 
conditions, fixing the value of w in Eq.Q, and searching 
for particular shapes of V n {x,t) = p n W(x ), n = 1,2,.. 
such that the stationary SE 

-dlS n /2p + p n W{x)S n =uS n , n = 0,1,2... (9) 
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has a solution S n (x,oo) which satisfies the boundary con¬ 
ditions 

S n [x,oo) ~ exp(±i^2poox), x —► ±oo. (10) 

The Sturmian eigenfunctions S n (also known as Stur- 
mians) differ for positive and negative w’s. As seen from 
Eq.(10), for oo < 0, all S n (x ) exponentially decay on both 
sides of the well, S n {x) ~ exp(— ^2p\oo\\x\) for |a;| —> oo, 
so that p n W(x) has a bound state at the chosen energy 
oo [ID]. For w > 0, a Sturmian contain outgoing trav¬ 
elling waves, S n [x,uo) ~ exp{±iy/2poo\x\) as x —> ±oo. 
This can only be the case if p n W(x ) is a complex valued 
emitting potential, which, in turn, requires Imp n > 0 for 
oo > 0. In general, as oo changes from —oo to +oo, a cho¬ 
sen p n {pj) traces a continuous trajectory in the complex 
p-plane. 

From Eq. © follows an orthogonality relation 

OS'mMISntu;)) = f S mM W(x)S mM dx =(11) 

fimn X (5 n (w)|5„(w)), 

where S mn is the Kroneker delta. The Sturmians are 
also known to form complete sets (we refer the reader 
to Ref.‘17) for a detailed discussion). Thus, to construct 
a physical solution T(x,t) describing particles which es¬ 
cape from the trap, we expand T(a;,a;) in © in the basis 
of S n , 


= y^B n (u])S n (x,u;), 


( 12 ) 


where the coefficients B n (oo) are to be determined. In¬ 
serting (121 into Eq. ([s]) , after adding and subtracting 
J2nPn WB nS ni we have 

+ (13) 

n k 

p n W(x)B n (u])S n (x,u])} = 0. 


In our quadratic case ©, multiplication of Eq.( |13[ ) by 
S m (x,oo) and integration over x, yields the following set 
of equations for B n (oo), 

M£l[v 2 B^ + (£-p m )B m }- (14) 

« 2 ^(2mWb; + mP|b b ) = o, 

n 

where a prime denotes differentiation with respect to oj, 
and 

M mn = J S m (x,uj)W(x)d J UJ S n (x,uj)dx = (15) 

(5 ro (u,)|5«M). 


Equations (14) and (PL5l) are the main achievement of the 


Sturmian approach: the problem of solving a partial dif¬ 
ferential equation ([2]) is reduced to solving a system of 
second-order ordinary differential equations. With only 
few terms in Eqs.(fl4|) usually needed, the Sturmian ap¬ 
proach offers a significant computational advantage in 


case of many dimensions HM3- In the one-dimensional 
case considered here, it can offer a further insight into 
the physics of scattering by time-dependent potentials 
and simplify calculations in certain limiting cases, as we 
will demonstrate next. 


IV. QUADRATIC ZERO-RANGE MODEL. THE 
ADIABATIC LIMIT 

We start with the simplest case of a zero-range poten¬ 
tial, 


W(x) = S(x), 


(16) 


which for £ — v 2 t 2 < 0 support a single adiabatic bound 
state ( 9(x ) = 1 for x > 0 and 0 otherwise) 

(j) 0 {x, t ) = [-2pE 0 (t)] 1/4 [6(x) exp(i^2pE 0 (t)x) + (17) 

0(-x) exp(-i^2pE 0 (t)x)\, 

with an energy 

E 0 (t) = -p{£-vH 2 ) 2 /2. (18) 

There is a only one Sturmian na, m 

Sq(x, uj) = [0(x) exp(i^/2p,ujx) + 0(— x) exp(—i^/2pcox)}, 

(!9) 

and the corresponding Sturmian eigenvalue, given by 


poM = i\J 2u/p, 


( 20 ) 


is single-valued on a two-sheet Riemann surface TZ of y/uj, 
cut along the positive semi-axis. With no other Sturmi¬ 
ans present, and = 0 since 5(0, w) = 1, taking 

complex conjugate of Eqs.([T4| [2D] yields 

B* 0 " + v~‘ 2 q 2 B* 0 = 0, q^)=^£-pi M (21) 


must be integrated along the contour running along the 
real w-axis above the cut the first sheet of 7 Z, where 
Sq(x,ui) satisfies the required outgoing/decaying waves 
boundary conditions (10). 


Equation (21), which is exact, can now be read in a com¬ 


pletely different manner. It has the form of a stationary 
SE describing a ” particle” of a ” mass” 1/2 with a ” coor¬ 
dinate” oj, of an ’’energy” £, scattered by a ’’potential” 
W(w) = Po(w), with v playing the of the ’’Planck con¬ 
stant” h ]21]. [We will always use the quotes when we 
refer to the fictitious ’’particle” in Eq. © , in order to 
distinguish it from the real particle described by the SE 
Q.] We note that the ’’potential” W, shown in Fig. 
2, has a valley (ReW < 0, ImW = 0) for oo < 0, and 
becomes purely absorbing, ( ReW = 0, ImW < 0) for 
oo > 0. 


Properties of equations of the type (21) are well known 
(see, e.g., [13]). As oo —► —oo, q(oo) —> oo, while W(oo) 
becomes flatter, W'(oo) ~ l/\/jw[ —> 0, so that B(oo) 
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can be expressed in semiclassical form ns, in terms of 
’’incoming” (+) and ’’outgoing” (-) ’’waves”, 


B* 0 M 


- 7 == exp[- f q{uj')dJ} + ( 22 ) 

\A/M vj 

A~ i 

— . exp]— / oo—oo 

L vj J 


where are unknown constants to be determined. We 
do not expect the particle to acquire very high energy, 
and must, therefore, require that B(ui) —> 0 as oo —> oo. 
Thus, taking the principal branch of the square root in 
Eq.© , we have 






exp[ 


|g(u; 7 ) | dco' / ] —> 0, oo —> oo. (23) 


Finally, inserting (12 1 and (22) into Eq.Q we note 
that as t —y ±oo the integral over to may be evaluated 
by the stationary phase method [23] , and as t —> — oo we 
have (details are given in the Appendix) 



CO 


4/(a;,t) ~ 2v\fn~iA + cj)(x,t) exp [—i 


E(t')dt'}. (24) 


This describes a particle trapped in the ascending bound 
state, which approaches the continuum threshold from 
below. Similarly, for t —> oo we have 


\k(x,t) « 2v\J—niA if>(x,t)exp[—i 


E(t')dt'\ (25) 


+<5'k(a;,f), 


where first term describes a particle trapped in the de¬ 
scending bound state, moving away from the continuum 
threshold. For the probability to complete the passage, 
and remain in the bound state we, therefore, have 


Tjstay 

M)0 


14-I 2 

i^+r 


(26) 


which is less or equal to one, as guaranteed by the ab¬ 
sorbing nature of the ’’potential” for x > 0 . 

Reduction of the original time dependent problem to the 
one of determining the reflection coefficient of a complex 
valued barrier allows us to prove the existence of the 
adiabatic limit in case the bound state ’’turns” without 
touching the continuum £ < 0. Now ’’absorption” repre¬ 
sents the loss of the particle to the continuum, and to be 
absorbed, the ’’particle” must first cross the ’’classically 
forbidden region” (see Fig. 2), impenetrable in the ’’clas¬ 
sical limit” v —> 0. This is the adiabatic theorem. The 
behaviour for £ > 0 requires somewhat more attention, 
and we will consider it next. 


V. A ZERO-RANGE WELL: JUST TOUCHING 
THE CONTINUUM 

With £ = 0, we have E 0 (t = 0) = 0, so the bound state 
of a zero-range well approaches the continuum thresh¬ 
old, and ’’touches” at the moment it ’’turns” to begin 


FIG. 2. (Color online). Schematic diagram for the zero-range 
model. A ’’particle” with an ’’energy” £ is scattered by a 
complex-valued ’’potential”, W, whose real and imaginary 
parts are shown by the solid and dashed lines, respectively. 
Also indicated is the ’’classically forbidden region”, separating 
a ’’particle” with £ < 0 from the ’’absorption” region w > 0 . 
The inset shows the energy of the bound state given by (18 1 , 
Eo(t) = p - 1 / 5 u - 4 / 5 .Eo(t), a s a function of. t = fj, 1 ^ 5 v 4 ^ 5 t. 


the downward leg of its journey. In this special case the 
equation for Bo, 

B'o + b^/ujB 0 = 0, b = \/2//i exp( 37 ri/ 2 )/ v 2 , (27) 


can be solved analytically in terms of the Bessel functions 
[241 . The solution which vanishes as ui —t oo, is given by 

B{w) = z = —e 3 ^ 4 ^ 4 , (28) 

where H„\z), j = 1,2 is the Hankel function of the j- 
th kind [25]. As u —t oo, we have (omitting inessential 
phase factors) 


B(uj) ~ u 1/8 exp [—(1 -I- i)Koj 5 / 4 ] u ^oo 0, (29) 


where K = 2 7 / 4 /5/i 1 / 4 > 0. Equation (29) is readily 
recognised as a special case of Eq. (23), with q(u) = 
(2 w/^i ) 1 / 4 exp( —in/ 4). To find the asymptotic form of 
B(co) for ui —» —oo, we use the formula connecting the 
values of Hy^ 5 (z) on the ray z' = pexp( 37 ri/ 4 ) with those 
along 2 = z' exp(* 7 r) = pexp(7ni/4) (see [2S], Sect. 3.62). 


H 2 /l( 2 ') = 2cos(27r/5 )H^l(z) + exp(-2ni/5)H^ 5 (z). (30) 


2/5 


2/5 ^ 


Recalling that h1 1,2 \z) ~ {2/nz) 1 ^ 2 exp [±*(2 — vn/2 — 
n/2)\ 55], and taking complex conjugate of Eq.(30), we 
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( 2 ) 

identify Hy 5 (z) with the ’’incoming wave” 
which gives 


in Eq.(22 I, 


P 0 S ‘“» = 4cos 2 (2tt/5) « 0.38197. (31) 


Thus, for a narrow well such that its bound state just 
’’touches” the continuum at t = 0 , the probability to 
remain in the well is independent of the rate of change 
of the potential. There is a perfect balance: a rapidly 
changing well is more likely to eject the particle into the 
continuum, yet the time the bound state spend near the 
threshold is short. If the well changes slowly, this time 
is longer, yet the particle is ejected less efficiently. As a 
result, there is no adiabatic limit as v —> 0 , and the value 
of Pg* ay is always given by Eq.(31). Below we will show 
that Eq. (31) has a more general meaning, also beyond 
the zero-range model considered in this Section. 


VI. A ZERO-RANGE WELL: THE GENERAL 
CASE 



Y 


FIG. 3. (Color online) The probability Pgg a!/ vs. 7 for the 
quadratic zero range model (161, obtained by integration of 
Eq.( |21| ) (solid), and by solving numerically the original SE 
Q (dots). 


No analytic solution of (27) is known (at least to us) for 


£ / 0, so the equation must be solved numerically. We 
note first that, for a narrow well (16), Pqq^ is determined 
by a single dimensionless parameter 


_ £p 2 / 5 

7 “ V 2 / 5 ' 


(32) 


Indeed, in the scaled variables r = and y = 

fj?/ 5 v 2 / 5 x, the SE pi reads id T ^(y, r) = — <9 2 'F/2-|- (7 — 
r 2 )<5(?/)'k, and Eq.([27|) only needs to be solved for v = 1, 
and various values of £ = 7 . The dependence of Pqq ( 7 ) 
on 7 is shown in Fig. 3. The probability tends to 1 for 
7 — > —00, where the ’’absorbing potential” in Fig. 2 
is separated by a broad ’’classically forbidden” region. 
At 7 = 0 the curve passes through the value given by 
Eq.(|3l|), P stoy (0) = 4 cos 2 (27t/ 5), and tends to zero as 
7 — > 00, i.e., when the ’’particle” can penetrate deep into 
the ’’absorbing region”, and nothing is ’’reflected”. 

Using Figure 3, it is easy to predict the behaviour of 
the retention probability P sta v as a function of v, for a 
given £. For £ < 0, and v « y\£ | 5 / 2 the passage will 
be adiabatic, with almost none of the particles lost. For 
£ > 0 and v « y\£\°J 2 , the bound state will disappear 
for a long time (see inset in Fig. 2), and none of the 
particles will be recovered when the it finally reappears. 
With v —> 00, 7 will vanish for any choice of £, and we 
have 



FIG. 4. (Color online)The probability Pq, g ay vs. v for the 
quadratic zero range model (|16[), for different values of £. 


the original SE ([ 2 ]). The results shown in Fig. 5 demon¬ 
strate that Po(t) = |(</>o(t)| v I , (t ))| 2 undergoes oscillations 
before reaching the asymptotic value Po(t) = Pgg ay , 
when the bound state is well removed from the contin¬ 
uum. 


VII. A RECTANGULAR WELL: THE 
ADIABATIC LIMIT 


lim v ^ 00 P^ ay {£,v) = 4cos 2 (27r/5) (33) 

so that a rapidly changing zero-range well will retain the 
particle in about 38% of all cases, regardless of the value 
of £. The dependence of (£,?;) on v for different 
values of £ is shown in Fig. 4. 

Finally, in order to study the evolution of the population 
P 0 (f) of the moving bound state, we solved numerically 


A more realistic case of a rectangular well of a width 

2a, 

W ( x ) = [6{x + a) 9 (a — x)\/2 a, (34) 

is somewhat more involved. There are two types of Stur- 
mians, symmetric and antisymmetric about the origin, 
S n {x, w) = S n {—x,ui), and T n (x,uj) = ~T n (-x,ui). For 
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FIG. 5. (Color online) The population of the moving bound 
state, -Po(f) = |(<Ml)|^(f)>| 2 , vs. T = p 1 ^ 5 v 4 ^ 5 t for three val¬ 
ues of 7 = £/r 2 / 5 v -2 / 5 . For 7 = 1 the bound state disappears 
at r = — 1 , and reappears again at r = 1 . 


— a < x < a , these are given by 

S n (x,u) = cos(p n x)/[I + sm(2p n a)/2p n } 1/2 , (35) 

n = 0,2,4..., 

and 

T n {x, uj) = sin(p„a:)/[l - sin(2p„a)/2p n ] 1/2 , (36) 

n= 1,3,5... 

so that (S n \S n ) = (T n \T n ) = 1. Since the matrix el¬ 
ements in Eqs.(14| couple only Sturmians of the same 
parity, (S m (w)|T^ (w)) = (T m (ui)\S^\w)) = 0, we may 
limit our analysis to the case where a particle is pre¬ 
pared initially in a bound state symmetric about the 
origin. The corresponding Sturmian eigenvalues p n , 
n = 0,2,4..., are then found by solving a transcenden¬ 
tal equation, 


sin (p n a )/cos (p n a) + ik/p n = 0, (37) 


where 


p n (co) = {2p[ui — p„(u;)]} 1,/2 , and k(u) = (2 pco) 1 ^ 2 . (38) 


Thus, p n (uS), is the magnitude of the rectangular poten¬ 
tial, real o complex, such that at a given energy w, there 
is a symmetric solution S n (x) of the SE ([9]), satisfying 
the boundary conditions (101. 

The interpretation of equations for 7?*(w) is similar to 
that given in Sect. IV. One may think of a fictitious 
’’particle” with an ’’energy” £ which can move on sev¬ 
eral complex valued ’’potential surfaces” >V„(w) = p* (w). 
On each ’’surface”, ’’absorption”, possible for uj > 0, 
accounts for the loss of the real particle to the contin¬ 
uum. There is also a possibility for hopping between the 


’’surfaces”, facilitated by matrix elements Mmn(u>) and 
If a particle is prepared in the rn-th state of 
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FIG. 6 . (Color online) Real (a) and imaginary (b) parts of 
the ’’channel potentials” = pa\V n vs. U = pa 2 uo for 
n — 0,1, 2. In a), also shown by a dashed line is the large-ta 
asymptote of Re Wo, oja — n 2 /2pa. 


the deep well, we must look for a solution of this ’’coupled 
channels problem” containing, as u> —> — oo, an incoming 
wave on the m-th ’’potential surface” and, possibly, ’’out¬ 
going waves” in all other ” channels”, 

A+ j 

B* n {u) ~ S mn — ^exp[- / q m {u)')dio']+ (39) 

v Qm V J 

A~ i 

n expf— / q n {uj')du> '1, w —>■ — oo, 

v J 

where m,n = 0,2,4.... The probabilities for a particle 
to start in the ?n-th, and end up in the n-th adiabatic 
bound states, are given by 

P ™ = \aB’ m,n = 0,2,4... (40) 

The ’’potentials” W’„(w) are shown in Fig. 6 for m = 
0,2,4. We have ImW n (±) = 0 for w < 0, where Stur¬ 
mians are just bound states of a real potential well. We 
also note that 

lim ReWniyj) = ua — (2 n + 1) 2 tt 2 /2pa, (41) 

to —^ioo 

and 


lim ImW n (uj) = 0. (42) 

uj— y oo 

This is because for a large negative u>, the Sturmians 
tend to the eigenstates of a potential box with infinite 
walls at x = ±a. Since the energy of the state is w, p n , is 
found by subtracting from u> the energy of the n-th state, 
as measured from the floor of the well. In the opposite 
limit, uj —^ oo, the particle becomes bound at the top 
of an infinitely high rectangular barrier. These bound 
states, quantised between the sharp potential drops at 





















VIII. A RECTANGULAR WELL: THE 
SINGLE-STURMIAN APPROXIMATION 
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FIG. 7. (Color online) Real (solid) and imaginary (dashed) 

(2) (2) _ o 

parts of the correction term Mq 0 = yaM 00 vs. lo = ya w. 


x = ±a, are essentially the same as those quantised be¬ 
tween the walls of an infinite potential box [5Hj . Since 
the Sturmians cease to depend on w as w -} ±oo, 

S n (x,Lo) ~ 9(x + a)6{a — x) cos[(2n + l)7r:r/2a], (43) 

the matrix elements, coupling the ’’potential surfaces”, 
vanish in the same limit, 

0, w -> ±00. (44) 

Both Mmn and Mmn are singular at the threshold w = 0, 
as is explained in the Appendix B. In particular, for 
( 2 ) 

Mq 0 , which will be required in the next Section, we find 
Mqq (ui) ~ w -1 ' 5 (see Fig. 7). 

We can now formulate the adiabatic limit for a particle 
prepared in the ground state of a rectangular well, m = 0, 
at t —► — oo, provided the state ’’turns” before reaching 
the continuum threshold, £ < 0. As in the case of the 
zero-range potential, the ’’absorbing region” is in Fig. 
6 separated by a ’’classically forbidden region” and be¬ 
comes inaccessible for a ’’particle” incident on the n = 0 
’’potential surface” as v —> 0. There is, however, a pos¬ 
sibility to access the ’’absorbing potential” in Fig. 6 by 
hopping to a different ’’potential surface”. But as v —> oo 
the hopping also becomes improbable, since the solutions 
on different ” surfaces” become highly oscillatory, and the 
integrals involving Mg n ’ (ui) vanish. We, therefore, have 
the adiabatic limit 

P^ n ay (£ <0,v^0)^5 On . (45) 

This result is easily extended to other initial states, m ^ 
0. The behaviour for other values of £ and v requires 
more attention, and we will consider it next. 


The above discussion suggests that if the trapping po¬ 
tential changes sufficiently slowly, one can largely neglect 
scattering into other bound states of the well, thus leav¬ 
ing a few, or indeed just one, equation in (14). For a 
particle arriving in the adiabatic ground state, m = 0, 
we, therefore, write 

[v- 2 (£-pl) + (S o \dlS o y}B* = 0, (46) 


B, 


where we have retained the diagonal correction term 
Mqq’Juj). With no analytical solution available for 
Eq.p6|), we have to solve it numerically. 

In the dimensionless variables r = 4 ma 2 t, y = x/2 a, 
[£ = 4 ya 2 £, and v = 8y 3 ^ 2 a 3 v the SE (p| reads 
id T ^{y,r) = -dfa/2 + \£ -v 2 T 2 )]6(y + 1/2)61(1/2 -y)T. 
and we must solve Eq.(W6|) for a particle of a unit mass 
in a well of a unit length, replacing to with ya 2 uj. The 
results for P^^are shown in Fig. 8 together with the ex¬ 
act curves, obtained by solving numerically the original 


SE (§. 

The exact results are worth a brief discussion. For for 
the ground state just ’’touching” the continuum thresh¬ 
old, £ = 0, the P c 
Eq.J3lb, 


QQ ay tends to the constant value in 


P ° l av (£ = 0,u -a 0) -a 4cos 2 (2tt/5). (47) 

This can be understood by scaling the variables in Eq. ([2]) 
in a different way, so as to put to unity the parti¬ 
cle’s mass y as well as v, i.e., r = i.e., 

y = y 3 / 5 v 2 / 5 x. With this we also have W ( y ) = [&(y + 
y 3 / 5 v 2 / 5 a)6(ii 3 / 5 v 2 / 5 a - y)\/2y 3 / 5 v 2 / 5 a. As v —► 0, the 
width of W(y) tends to zero, and we recover the zero- 
range result ( |47| , which holds universally for all values 
of v and, in particular, for v = 1 . 

For a rapidly changing rectangular trap, v —> oo, the 
particle always returns to the well, regardless of whether 
the adiabatic state ’’turns” before touching the contin¬ 
uum, just touches it, or even disappears for a while. An 
yet different type of scaling can be used to explain why. 
Putting to unity the well’s width as well as v, we have 
a particle of a mass ~p = 4 a 2 v 2 ^ 3 y, and a new param¬ 
eter, £ = £/v 2 / 3 . As v —» oo, we have a picture of a 
very heavy particle, y —> oo, brought to the continuum 
threshold, £—>(). and then down again. The massive 
particle has no chance to escape, and we have 


Poo av (£, v -A oo) -A 1, (48) 


which holds for all finite values of £. 

Finally, if a bound state disappears, the particle’s state 
is a wave packet of continuos states, which spends a du¬ 
ration of 2^/~£/v spreading away from the region. For 
v —> 0 the time of spreading is very long, so that little is 
recaptured after the bound state reappears at t = y/£/v. 
Thus we have 


Poo av {£ > 0,w -A 0) -a 0. 


(49) 











FIG. 8. (Color online) Probability to remain in the ground 
state of a rectangular well, P 0 S q“ s vs. v = 8 fi 3 ^ 2 a 3 v for 
£ = —4, — 2,0, 2, 4 (solid). Also shown are the single-Sturmian 
approximations to these probabilities obtained with (dashed) 
and without (dot-dashed) the last term in Eq. (461. The in¬ 
set shows the total probability Pg tay (thick solid), the exact 
(solid), and the single-Sturmian approximation to Pq^ 
(dashed). The vertical dashed lines indicate the three regimes 
described in Sec. VIII. 


pstay 

M) 


The single-Sturmian approximation for Pqq 0!/ , obtained 
by solving Eq.(46), is in good agreement with the exact 
result for v < 40. Comparing the two curves with the 
total probability to stay in the well, P^ tay = Po^ v , 
shown in the inset in Fig. 8 helps identify three approx¬ 
imate regimes. 

A) Slow passage. For v < 5, we have P^ tay « P^ ay . 
The loss and recapture of particles is determined by inter¬ 
action of a single bound state with the continuum. There 
is no scattering into other bound states. Mathematically, 
the problem reduces to solving a single equation (46) (see 
Fig. 9a). 

B) Intermediate passage. For 5 < v < 80, we have 
> P°l av , ( P° tav - P^ ay yP stay « 1, with P^ ay 

is correctly described by Eq.(46|. This suggests that a 
downward bound initial state recaptures some of the par¬ 
ticles, and later each new bound state, which enters the 
deepening well, scoops some more (see Fig. 9b). This 
regime can be described by solving Eqs.([T4| iteratively, 
using the solution of (46) as an initial approximation. 

C) Rapid passage. For 80 v we find notable dis¬ 
crepancies between the single-Sturmian approximation 
for Poo Qy , and the exact result. This indicates that the 
loss to continuum is accompanied also by transitions be¬ 
tween different bound states. Mathematically, this re¬ 
quires solution of the full ” coupled channels problem” 
(14) (see Fig. 9c). We note that in the case of several 
spatial dimensions reduction of the original problem © 
to that of solving a system of ordinary differential equa¬ 



FIG. 9. (Color online) Rectangular barrier: p° tay V s. r = 
4 ga 2 t for v = 0.05; b) same as a), but for v = 10. Vertical 
dashed lines indicate the moments when the first and the 
second excited states enter the well: c)same as a), but for 
v = 200. 


tions may be a significant simplification. 


IX. UNIVERSALITY OF THE ”38% RULE” IN 
THE v -> 0 LIMIT 


We have shown that in the two cases considered above, 
there is no conventional adiabatic limit for a ground state 
just touching the continuum threshold. Rather, as v —> 0, 


the probability to remain in the state is given by Eq.(311, 


and equals approximately 38%. It easy to show that, 
for a quadratic evolution ([6]), this result holds true in 
one dimension for a finite-width potential of an arbitrary 
form. Indeed, scaling the time and coordinate so as to put 
to unity v and the particle’s mass p 1 while maintaining 
the normalisation f W(x)dx = 1, 


r = /WV/V, 


y = M 2/ V /5 x, 


(50) 


converts the SE © (lA°) = 0) into 

id T ^(x, t) = -d 2 ^/2 + T 2 W{y)' L, (51) 

where W[y) = y~ 3 ^v~ 2 ^W{p~ 3 ^v~ 2 ^y). As v —> 0, 
we have W(y) S(y) for any choice of W{x), so that 
Poo“ y is given by Eq. (31). To illustrate this, we plot¬ 
ted the results for a rectangular well (34), and a cut-off 
parabolic potential 


W{x) = 1.5 [6{x + a) 6 (a — x)\x 2 /a 3 , 


(52) 


in Fig. 10a. 

The case where the m-th excited state of the well touches 
the continuum requires more attention. Let the evolution 
of the potential be such that at t = 0, in the potential 
-V<%( x), we have E m — 0. Returning to Eq.(|14| we 
note that as v —> 0, the last sum in it may be neglected. 
Also, in this limit ’’absorption” of the ’’particle” occurs 
in a small vicinity of w = 0. Thus, if we can demonstrate 
that, for co « 0, B m (uj) satisfies Eq.(27), obtained ear¬ 
lier for a zero-range well, we will have also proven that 
Pmm( v -A 0) = 4cos 2 (27r/5), for all m’s and all potential 
shapes. 
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where C > 0 is a real constant, and E m (p°) = 0. Solving 
the pole condition ik + n(p) = 0 for p , we have 


Pm(u) = i\/2puj/C m + p ° m , 


(55) 


where we recalled that our derivation is for a particle 
prepared in the m-th state of the deep well, and added 
the index to, where required. Inserting (55) into Eqs.(141, 
neglecting all but one of them, and noting that = 
Pm, for u) « 0 we have 


B„ 


+ cyfu}B m = 0, c=-iy/2p,/c m v 2 (56) 


The similarity between Eqs. (27) and (|57j) in the region 
of interest allow us to conclude that for any state <j) m , 
m = 0,1, 2.., ’’touching the continuum threshold” 


lim v ^ 0 P^y = 4cos 2 (27r/5)5 mn »= 0.3819 x <5 mn .(57) 


FIG. 10. (Color online) a) the probabilities to remain in 
the ground state of the square (dashed) and parabolic (solid) 
wells, Pq o° v for £ = 0, vs. v = 8p 3 ^ 2 a 3 v The inset shows dif¬ 
ferent shapes of the parabolic well ( |52| ); b) same as a), but for 
the first 4 states of a rectangular well ’’touching the contin¬ 
uum threshold”. A dashed line shows P 0 S Q ay for £ = —0.015. 


This general result is valid for any potential, provided 
the scattering length C has a simple pole when the bound 
state <p m joins the continuum, C(p) ~ 1 /(p — p°), where 
E m (p°) = 0. This condition is fulfilled, for example, for a 
rectangular well (34), with the results for various excited 
states, obtained by integrating Eq. {46), shown in Fig. 
10b. 


To show that this is the case, we use the standard ap¬ 
proach, commonly used to describe potential scattering 
at low energies 22]. To this end, we will consider a 
potential pW(x), whose m-th bound state lies just be¬ 
low the threshold, construct a continuum state for a 
small positive energy, u > 0, and look for the con¬ 
dition under which the scattering amplitude diverges, 
S(E,p ) —> oo. Namely, our state will have the form 
exp (±ikx) + S(ui, p) exp(^fikx), k = yJ2pui, for x>a. In 
the low energy limit, ka << 1, the wave length of the 
particle is large, and the well is characterised by a sin¬ 
gle parameter, the logarithmic derivative of the bound 
state’s wave function at x = ±a, which we denote —k, so 
that (p m (x,p) ~ exp[— n{p)x\. We note that k depends 
on the potential shape via p, but not on the energy w, 
as long as u> is small. Matching the log-derivatives at 
x = ±o, and using ka » 1 then yields 


S{E,p) 


ik — n(p) 
ik + n(p) ’ 


(53) 


which diverges whenever ik + n{p) = 0. The condition is 
usually used to obtain the pole in the complex cu-plane, 
given a real value of p [22] . We, on the other hand, require 
the value of p, given a real value of ui, and need to make 
an additional assumption about how n depends on p. 
The scattering length, defined as C = — 1/k, is known to 
remain real, diverge, and change its sign as the shallow 
bound state moves toward the continuum threshold, and 
eventually becomes a virtual state [52]. Thus, we will 
assume k to be a linear function of p, 


«(p) ~ C(p- p°), (54) 


X. SUMMARY AND CONCLUSIONS. 


In summary, we have analysed, in one dimension, the 
evolution of a particle prepared in a bound state of a 
trapping potential, whose magnitude has a simple maxi¬ 
mum at t = 0, as described by Eq. (§• There are three 
possible scenarios for the state, which first approaches 
the continuum threshold, and then moves away from it. 
It may (i) turn before reaching the continuum threshold, 
(ii) just touch it once, or (iii) cross the threshold and 
temporarily disappear. Whether the particle remains in 
the trap, or is lost to the continuum, depends on how 
fast is the variation of the trapping potential. 

In the slow passage limit, the particle always remains in 
its initial (mth) state, provided the state ’’turns” before 
reaching the threshold, in accordance with the Adiabatic 
Theorem. If the state touches the threshold, the proba¬ 
bility to remain in it Pmm is approximately 38%. This 
result holds universally for all excited states and various 
potentials, under a very general assumption about the 
behaviour of the scattering length (541, and replaces the 


conventional adiabatic limit. If the bound state disap¬ 
pears for while, a particle ejected into the continuum has 
sufficient time to move away from the potential. Thus, 
there is 100% loss to the continuum, and nothing is re¬ 
covered when the state reappears. 

In the rapid passage limit, the outcome depends on the 
choice of the potential. Thus, for a zero-range well, Pqq ay 
tends to the same 38% limit, regardless of whether the 
bound state turns, touches the threshold, or crosses it. 
This appears to be a consequence of a perfect balance 
between the time a bound state of a (5-well spends near 
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the threshold, and the efficiency with which the particle 
is ejected. On the other hand, in the case of a rectangular 
potential, a rapidly evolving well always retains the par¬ 
ticle in its original state, whichever the fate of the bound 
state. 

The general case of a passage which is neither slow 
nor fast is conveniently studied in the Sturmian repre¬ 
sentation. Unless the potential changes very rapidly, it 
is sufficient to employ only one Sturmian state, and the 
task of solving the time-dependent SE ([2]) reduces to that 
of evaluating the reflection coefficient of a complex valued 
’’potential”, where absorption of a fictitious ’’particle” ac¬ 
counts for the loss of the real particle to the continuum. 
For larger values of v, several Sturmian states need to be 
taken into account, and the picture is that of a ’’parti¬ 
cle” capable to moving on several absorbing ’’potential 
surfaces”. In general, one can loosely identify three dif¬ 
ferent regimes. If the passage is sufficiently slow, the 
state ejects the particles on its way up, and then recov¬ 
ers some of them on its way down. For faster variations, 
the original state recovers its share of the particles, while 
more particles are scooped by other states, which enter 
the well as its depth increases. At yet larger u’s, the 
loss to the continuum is accompanied by scattering into 
other bound states, and one needs to solve a full ’’coupled 
channels problem” (141. 

Verification of the above theory is within the capa¬ 
bilities of modern experimental techniques, e.g., of the 
laser-based methods for containing cold atoms in quasi- 
one-dimensional traps. In spite of a practical difficulty of 
assuring that the state just touches the threshold, this re¬ 
sult should be amenable to an experimental verification. 
Figure 10b shows PgQ ay for a state that turns shortly be¬ 
fore reaching the continuum, 6 = —0.015, closely follows 
the 8 = 0 curve before shooting up to its adiabatic limit 
P o r = 1 for very small values of v. Thus, the condition 
8 can be fulfilled approximately, provided v is chosen to 
be not too small. 

Among other advantages offered by the Sturmian tech¬ 
nique is a simple interpretation of the adiabatic condition 
for a state which turns before reaching the threshold. In 
this case, in order to be absorbed the fictitious ’’particle” 
must first cross a ’’classically forbidden region” in Fig. 
2. With v playing the role of a ’’Planck’s constant”, this 
becomes improbable, if the passage is slow. How P^m 
tends to the adiabatic limit as v —>• 0 can then be studied 
by evaluating the corresponding phase integrals. We will 
consider this in our future work, together with extending 
the analysis to several spatial dimensions, different tem¬ 
poral evolutions, and the case of several identical bosons 
trapped in the same bound state. 
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XII. APPENDIX A. 


For a t > 0, the stationary phase approximation to the 
integral ([7]) evaluated along the contour specified in Sect. 
IV is given by 

/•UJ 

u>, x) exp[— iuit + i / q(u>')du)']du> 

J r J ujq 

~ [27r*/g(w s )4'"(w s )] 1/2 S'(w s ,a;)exp[*$(a) s ,t)], (58) 

where $(w, t) = — iu>t + i q(ui')doj' , and u s < 0 is 
defined by 


I(t) EE / q (u)- 1/2 S( 


q(v s ) = t. (59) 

Given the time evolution of the magnitude of the 8- 
potential, there are three quantities, each of which can 
be used as an independent variable. These are the time 
itself, r, the well’s depth V(t) = So — v 2 r 2 , and the 
energy of the adiabatic bound state supported by the 
well, E(t) = — /W 2 (t)/ 2. It is readily seen that q(w) in 
Eq.([2T|) gives the time r, at which E[t(w)) = w 

q(u>) = t(u>) = [8 — / qi\ 1 / 2 /v (60) 


Let the lower limit in the integral in the exponent in 
Eq.(58) be uj 0 = — g,8 2 /2 if 8 > 0, and 0 otherwise. This 


ensures that q{u>) is always real non-negative for u) < 0. 

Changing variables w-► r(w), and integrating by parts, 

we have 

[ Tduj/drdT = tE(t)\ t t Z t t {^ o) - [ E(r)dr. (61) 


With t(ui s ) = t, and either t(u>q) or E(t(u!q)) vanishing, 
we have 


<h(w s ) = - f E(t )dr. (62) 

JlJq 

For the second derivative of the phase, $"(w s ), and the 
pre-exponential factor, we obtain 

<I>"(w s ) = q'(u s ) = [d£ , (r)/dr| r=t ] _1 = (63) 

— [2/j,v 2 t(fl 0 + v 2 t 2 )]~ l = — [2v 2 t\/— 2/j,E(t)]~ 1 . 
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g(uj s ) = S(x,u> s )/v / q(u> s ) = S(x,E(t))/y/t (64) 


Inserting (62), (63) and (64) into ( [58] ) , and using ( fl7| ) 
yields the term which multiplies A + in Eq. (25). Equa¬ 


tion (24) for a t < 0 can now be obtained as complex 


conjugate of Eq.(60). 













11 


XIII. APPENDIX B. 


For a rectangular potential of a unit width, a = 1/2, 
and a particle of a unit mass, /i = 1, we have 


(S m (w)|S n (u/)) = (65) 

F({Pm + Pn)/ 2) + F((p m - p n )/2) _ ^ x 

[F(p m )F(p n ))^ - 


where p m (u) = [2(u - p m (w)] 1/2 , p n {u') = [2(u - 
/^(o/)] 1 / 2 , and F(x) = Thus, the coupling ma¬ 

trix elements are given by 


(1) _ dG(p 

mi Pn) dp n 

mn ~ d Pn du' 


( 66 ) 


and 


„'=« ( 67 ) 


,,(2) = 9G{p m ,p n ) d 2 p n d 2 G(p m ,p n ) / dp n 

mn d Pn du' 2 dp 2 n \du' 

The divergencies of Mhh ’ at u = 2 come from the deriva¬ 
tives of p n , which has a branching singularity at u' = 0. 
It follows from Eq.(37) that, as u' —> 0, po ~ w ,1//4 , and 
Pn^o ~ o/ 1/2 . Therefore, for w->0we obtain [Mmh = 0 
since (S m (u)\S n (u)) = 1] 

M mo ~ w' 3/4 , ~ W -V 2 / or m ^ o, n. (68) 


Similarly, since the first term in Eq. (67) vanishes for n = 

TO, 


M, 


( 2 ) 

00 


' u 


- 1.5 


and ~ u 1 for to ^ 0. (69) 


For to ^ n the first term in Eq. (67) dominates, which 
leads to 


M md ~ w 175 > and M&1 ~ u 15 /or to ^ 0. (70) 
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